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We formulate a model for magnetic and superconducting ordering at LaAlOa/SrTiOs interfaces containing 
both localized magnetic moments and itinerant electrons. Though these both originate in Ti 3d orbitals, the for- 
mer may be due to electrons more tightly-bound to the interface while the latter are extended over several layers. 
Only the latter contribute significantly to metallic conduction and superconductivity. In our model, the inter- 
play between the two types of electrons, which is argued to be ferromagnetic, combined with strong spin-orbit 
coupling of the itinerant electrons, leads to magnetic ordering. Furthermore, we propose a model for interfacial 
superconductivity, consisting of random superconducting grains in the bulk STO driven, via coupling to the 
interface conduction band, towards long-ranged or quasi-long-ranged order. Most interestingly, the magnetic 
order and strong spin orbit coupling can lead in this manner to unconventional interfacial superconductivity, 
yielding a possible realization of Majorana physics. 



I. INTRODUCTION 

It was recently discovered that, although LaAlOa and 
SrTiC>3 are both insulators, the interface between them 
is metallic . Furthermore, the electrons at this interface 
have shown a variety of remarkable properties, including 
magnetism ~ and superconductivity . Magnetism and super- 
conductivity often appear in the phase diagram of strongly- 
correlated materials, where they compete. However, it is very 
unusual for them to occur simultaneously, which appears to 
be the case at the LAO/STO interface 4 '^' 7 . A basic question, 
then, is whether the same electrons are exhibiting both super- 
conductivity and magnetism or if, instead, there is a precise 
sense in which there are two species of electrons - two differ- 
ent electronic bands, for instance - one of which is supercon- 
ducting and the other of which is magnetic. In the former case, 
the superconductivity must be exotic, perhaps p-wave super- 
conductivity or a Fulde-Ferrel-Larkin-Ovchinnikov (FFLO) 
state 8 ' 4 . This would contradict the conventional wisdom that 
the interface electrons are simply exhibiting the superconduc- 
tivity of doped SrTiOs, which is presumed to be a phonon- 
mediated s-wave superconductor. 

The experiments that reveal interesting magnetic behavior 
fall into two classes, those which deduce magnetism from 
transport and those which attempt to measure it more directly. 
The former include experiments that observe hysteresis in the 
electrical resistance as a function of magnetic field ■ , which 
show that there is a magnetic field-driven first-order phase 
transition which has a large effect on the resistance. A mag- 
netic transition is the most natural hypothesis. These signa- 
tures are found up to temperatures in excess of 200K. The lat- 
ter include a torque magnetometry measurement that shows 
that a field as small as a few milliTesla leads to a large mag- 
netization, approximately 0.3 /is per interface unit cell. This 
implies that the system has ferromagnetic domains which be- 
come aligned by even a very small field. This experiment 
shows that the magnetic moment, which points in the plane, 
has an onset temperature that is at least as high as 40K and per- 
sists below the superconducting T c . Finally, scanning SQUID 



magnetometry finds micron sized ferromagnetic domains in 
a paramagnetic background. From their estimates, most of the 
interfacial electrons which are predicted by polar catastrophe 
arguments are paramagnetic. An order of magnitude smaller 
number of electrons are in ferromagnetic regions, and a two 
orders of magnitude smaller number of electrons are in su- 
perconducting regions. There does not appear to be any cor- 
relation between the magnetic and superconducting regions 
(unless there is spatial segregation in the z direction), which 
implies that magnetism in a region does not prevent super- 
conductivity from occurring in that region, counter to conven- 
tional wisdom. Taken together, these experiments constitute 
strong evidence that there are ferromagnetic domains at the 
LAO/STO interface that strongly affect normal state transport 
and are also present in the superconducting state. 

Since neither LAO nor STO is magnetic, there is clearly a 
puzzle here: what is the cause of (at least local) ferromag- 
netism at their interface? It has been suggested that there 
is a narrow band at the interface which gives rise to itin- 
erant electron ferromagnetism. Alternatively, the magnetism 
may be due to localized electrons which don't participate in 
the metallic (or superconducting) behavior. To make matters 
even less clear, there is evidence 1 for strong spin-orbit cou- 
pling due to the broken inversion symmetry of the interface 
(Rashba spin-orbit coupling), which would ordinarily be anti- 
thetical to a uniform ferromagnetic moment. 

An equally vexing problem is how ferromagnetism can co- 
exist with superconductivity 14 . Even if there were a sense in 
which different electronic bands were becoming ferromag- 
netic and superconducting, one would expect the magnetic 
moments of the former electrons to destroy superconductiv- 
ity in the latter. One possibility is that the system is in an in- 
homogeneous superconducting state, e.g. the FFLO state, as 
suggested in Ref. . Or, one could imagine larger scale inho- 
mogeneity, so that the system breaks up into domains, some 
of which are superconducting while others are ferromagnetic. 
But the ferromagnetic moment would be anti-correlated with 
superconductivity in either type of inhomogeneous state. For 
instance, in the FFLO state, the system forms superconduct- 
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ing stripes, separated by magnetic ones. This disagrees with 
the experimental finding of Ref. . Therefore, it is natu- 
ral to consider, instead, a p-wave superconducting state, but 
this begs question of what the superconducting mechanism 
is. Presumably, p-wave superconductivity must be due to 
electron-electron interactions, rather than the electron-phonon 
coupling. Thus, superconductivity at the LAO/STO interface 
is a puzzle. It is generally assumed that it is related to su- 
perconductivity in doped STO, but this assumption does not 
provide any clues to how it can coexist with ferromagnetism. 

In this paper, we present a physical picture for magnetism 
and superconductivity at the LAO/STO interface. Density 
functional theory calculations show that there are t2 g bands 
at the interface, corresponding to Ti 3d orbitals, with d xy and 
d X z,yz symmetry respectively ■ . We hypothesize that the 
d xy band forms a band of localized electrons which accounts 
for most of the charge required by the 'polar catastrophe' . 
According to our picture, the d xz , yz bands form metallic 
bands of itinerant electrons at the interface. Coulomb inter- 
actions between localized and itinerant electrons generates a 
ferromagnetic interaction between them, thereby leading to 
a ferromagnetic Kondo model - but one in which the itiner- 
ant electrons have significant Rashba spin-orbit coupling. By 
analyzing a spin-orbit-coupled ferromagnetic Kondo lattice, 
we argue that the localized electrons develop magnetic order. 
Thus, in our picture, the magnetic moment of the system is 
due primarily to localized electrons. 

Our picture for superconductivity is the following. We sup- 
pose that there are droplets of local superconductivity in the 
STO substrate. If the STO were doped, then these droplets 
would grow and percolate across the system, giving rise to 
superconductivity. In the absence of doping, this cannot oc- 
cur, and the STO substrate is insulating. However, these 
droplets can interact with the itinerant electrons at the inter- 
face. Through the proximity effect, superconducting droplets 
in the STO substrate can induce a gap in the itinerant elec- 
trons at the interface. Thus seeded, the itinerant electrons at 
the interface can develop long-ranged or quasi-long-ranged 
superconducting order. However, these itinerant electrons 
must move in the magnetic background created by the local- 
ized electrons. Naively, the magnetism should destroy the 
superconductivity. However, the strong spin-orbit coupling 
of the interface electrons allows these two competing phe- 
nomena to coexist peacefully . Spin-orbit coupling mixes 
s-wave and p-wave superconductivity, so that the s-wave su- 
perconductivity which is present in the droplets in STO can 
induce a mixture of s-wave and p-wave superconductivity 
at the interface. This mixture can tolerate a magnetic mo- 
ment, unlike pure s-wave superconductivity. This is very 
similar to the situation in proposals of topological insulators 
in contact with s-wave superconductors ~ , superconductor- 
semiconductor heterostructures • , and spin-orbit-coupled 
quantum wires 1 '-"'. Therefore, according to our theory, even 
though superconductivity is s-wave in STO, interfacial super- 
conductivity is unconventional, as a result of magnetism and 
spin-orbit coupling. 

An especially exciting consequence of the unconventional 
nature of the interfacial superconductivity is the possibility of 



realizing Majorana fermion physics. Indeed, the fabrication of 
narrow quasi-one dimensional conducting channels on an oth- 
erwise insulating LAO/STO interface 26 " 28 by 'writing' them 
with an atomic-force microscope (AFM) tip suggests a natu- 
ral implementation for the proposals of Refs. and 25. In 
the following sections we show that our model naturally gen- 
erates a 'helical' interfacial band structure that can be driven 
into a topological phase when in proximity to ordinary s-wave 
superconductivity. As shown by Kitaev [29], ID spinless su- 
perconductor with p-wave (or effectively p-wave) pairing sup- 
ports Majorana zero-energy modes at the ends. Although our 
droplet model can at best lead to quasi-long ranged supercon- 
ducting order, signatures of Majorana physics still remain, as 
shown in Refs. [30-32]. 

In Section II, we set up and justify a model of spin-orbit- 
coupled itinerant electrons interacting with Kondo spins. In 
Section 111, we give a saddle-point analysis of a large- N limit 
of this model. In Section IV, we solve this model numerically 
by the density-matrix renormalization group (DMRG) in the 
one-dimensional limit. In Section V, we analyze the super- 
conducting proximity effect in the presence of a helical wire 
due to the presence of ferromagnetic Kondo interactions and 
spin-orbit coupling. We show that such an interfacial super- 
conducting state might support Majorana zero-energy modes. 
We conclude in Sec. VI with the discussion of our results and 
proposing a schematic phase diagram for LAO/STO interface. 



II. SPIN-ORBIT-COUPLED FERROMAGNETIC KONDO 
MODEL 

We will now argue for a model of spin-orbit coupled 
electrons interacting with localized spins to describe the 
LAO/STO interface. The emerging consensus is that the elec- 
trons active at the LAO/STO interface come from the t^g 
bands of Ti 3d orbitals in the STO 3 ' 10 ' 16 ' 33 ~ 3S . Furthermore, 
band structure calculations and density functional theory 
suggest a picture of successive Ti 3d sub-bands near the in- 
terface being occupied as LAO thickness is increased (or gate 
voltage decreased). 

Due to the inversion symmetry breaking at the interface, 
the lowest sub-band is predicted to be d xy . It contains most of 
the charge required by the polarization catastrophe, but these 
electrons are thought to be localized, as seen from the low mo- 
bile carrier densities extracted from Hall transport measure- 
ments. We thus model these electrons as localized spins. Fur- 
ther indirect evidence for a picture involving localized spins 
in LAO/STO comes from transport experiments on a related 
system consisting of pure STO, with the doping effect of the 
LAO simulated by an polarized gel overlayer . The longi- 
tudonal resistance in this system exhibits a Kondo minimum 
as a function of temperature, indicative of the presence of lo- 
calized impurity spins. 

The lowest extended interface states are thought to have d xz 
and d yz symmetry. Due to the geometry of their orbitals, their 
Fermi surfaces are highly anisotropic, with heavy and light 
carrier directions. Spin-orbit coupling also plays an important 
role in determining the electronic band-structure, with the au- 
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thors of Ref. arguing for an atomic spin orbit (ASO) effect 
of about lOmeV. Furthermore, there is a Rashba contribution 
arising from the broken inversion symmetry of the interface 
and gating. Its magnitude is expected to be dependent on the 
gate voltage and the details of the sample, but in Ref. 11a 
value of a s» 10 — 50mcVA is obtained through a fit to a 
weak anti-localization measurement. (Larger values close to 
50meVA are obtained for larger gate voltages; smaller values, 
for smaller gate voltages.) The Rashba nature of the coupling 
was deduced from the dependence of the spin relaxation time 
on the elastic scattering time. The electronic band-structure 
for the two dimensional interface depends on the precise ratio 
of Rashba and ASO coupling. 

Motivated by recent experiments fabricating one dimen- 
sional conducting channels on otherwise insulating LAO/STO 
interfaces 27 ' 28 , as well as by our desire to realize Majo- 
rana physics, we find it useful to examine one-dimensional 
channels at the LAO/STO interface. In a ID channel, the 
anisotropy of the d xz , d yz orbitals suggest that for a very nar- 
row channel along the x direction, transport should primarily 
be through the d xz states . Therefore, the neglect of the d yz 
band is justified in this situation. In Ref. 27 and 28, con- 
ducting channels of thickness ~ 10 nm are constructed. For 
an effective mass of order the electron mass, such channels 
correspond to a transverse confining energy of ~ 2 meV, and 
hence in principle we need to consider many subbands at the 
energy scales in which we are interested. However, if nar- 
rower conducting wires could be constructed (e.g. 2 nm, cor- 
responding to ~ 40 meV), a one dimensional model would be 
more readily applicable. 

We are thus led to a analyze a minimal model of a single 
spinful band interacting with a large density of localized spins. 
The localized spins, which give the dominant contribution to 
the magnetic moment come from the d xy band. However, 
their tendency to order is driven primarily by their interaction 
with d xz , d yz electrons. Although the experiment of Ref. 39 
suggests an antiferromagnetic coupling between conduction 
and impurity spins, we believe that it is more natural, accord- 
ing to Hund's rule, take the interaction between d xy electrons 
and d xz , d yz electrons to be ferromagnetic. In either case, the 
Hamiltonian takes the form 



k x ,A,A ; 

+ J Y cL+a.A^A.A'Cfc^A' • S(q) 

k x ,q,\,\' 



a) 



b) 



(i) 



where cr, is the Pauli matrix acting on a spin degree of free- 
dom, S corresponds to classical spins due to localized mag- 
netic impurities, and we have suppressed spin indices in the 
notation. In a lattice version of this Hamiltonian the kinetic 
term is modeled by a nearest neighbor hopping t, whose mag- 
nitude is determined by band structure calculations ■ to be 
about 0.2 eV. From magnetoconductance experiments"' 12 we 
take a s» 10 — 50meVA, which translates to a lattice spin 
orbit coupling ctaT 1 ~ t/20 or even as large as t/4 (see be- 
low). The final term is the interaction between the itinerant 
electrons and localized spins. Although we believe that it is 





FIG. 1. (a) Sample one dimensional electronic band-structure. We 
expect a magnetic instability at ordering wave-vector q = 0, where a 
gap can open in the spectrum at the crossing point between the two 
bands. When this occurs, the band structure takes the form in (b) 



more natural, according to Hund's rule, to have a ferromag- 
netic coupling, J < 0, we will also, for completeness, con- 
sider the case of antiferromagnetic coupling, J > 0, with a 
magnitude of roughly one third of the bandwidth, as obtained 
in Ref. 39. 

Motivated by recent SQUID and torque magnetometry 
experiments, which find at least local ferromagnetism, we 
are interested in examining the magnetic instabilities of the 
Hamiltonian (2). The one dimensional spin-orbit coupled 
electronic band-structure in figure 1 shows that it is natural 
to expect ordering at wave-vector q = 0, since it is there that 
a gap can be opened and the electronic energy lowered. Of 
course, ordering could also occur at the wavevectors that con- 
nect various pairs of Fermi points. In the next section we will 
study a large-A^ limit of this model and find a tendency to- 
ward in-plane ferromagnetism, consistent with [ and 5]. In 
the following section we will reach a similar conclusion with 
DMRG. 



III. LARGE-TV ANALYSIS OF THE MODEL 

We now analyze magnetic instabilities of a generalization 
of the Hamiltonian (2) to N species of fermions c£ , with 
a = 1,2,..., TAT: 

H = Y [ e ( k <") C ktx C k x ,X S X,X> + ak x a X,X> C kl ,A C L,A' 
k x ,A,A ; ,a 

+ J Y C kl+ q ,\ S *,v4 x ,\> ■ Sq 
k x j^jAjA' 

= Y c £xM(k x ,S) xx ,clx (2) 

k x ,\.\' ,a 

where the matrix M(k x , S) is given by 

M(k x ,S)-- 



s(k x ) + JS z iak x + J(S X + iS y ) 

iak x + J(S X - iS y ) e(k x )-JS z 

(3) 
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We now integrate out fermions in Eq.(2): 

exp (-S cS [S}) = f DcdcU- s °&-f dT ^< ia ^-V 

_ e -Sn[S f ]-Artrln(iw„-M(fc x ,5)) 

(4) 

where Sq[S] is the Berry phase term for the localized spins. 
In the large-TV limit, the functional integral 



Z[S] = DS cxp (-S [S) - Mr In iw n ~ M{k x ,S) 



(5) 



is equal to its saddle-point value. The saddle-point equations 
are given by: 



d_ 
dS 



(s [s) 



Mr In 



iu n - M(k x , S) 



0. 



(6) 



The Sq[S] term is O(l) and is much smaller than the sec- 
ond term which is O(N), so it can be neglected. Since the 
temperature of interest is much lower than the energy scales 
associated with the couplings in (2), we can effectively set it 
to zero and convert Matsubara sum to an integral. Thus, we 
arrive at the following mean-field equations: 



d f du dk x 



dS J 2tt 2tt 



log [( iu -E+(k x ))(icu- E- (k x ,S) )] = 0, 



where E ± (k x ) are the two eigenvalues of M(k x , S): 
E ± (k x ,S)=s(k x )±Q(k x ,S) 



Q{k x ,S) = yf J* SI + J 2 Sl + (ak x + JS y Y 
The explicit evaluation of derivatives yields 



(7) 

(8) 
(9) 



/ 



dui dk x dQ(k x , S) 
2ir 2tt OS 



1 1 

icj-E+(k x ,S) iuj-E-{k x ,S) 

(10) 



The S'-dependent contribution to the energy of a particular 
background spin configuration at one loop reads 

/dk 
- (e+q(-e+) + £ fe -e(-£ fe -)) (13) 

To diagonalize, we take nearest neighbor hopping with ampli- 
tude normalized to 1, resulting in a kinetic term 



= 0. 



One can notice that we only obtain a non-zero contribution 
if E + (k x , S) and E~(k x , S) have opposing signs, so that the 
saddle point equation reduces to 



~27 



e(-E + (k x ,S))-@(-E-(k x ,S)) 



dQ(k Xl S) 



OS 



0. 



(11) 



Here Q(x) is the unit step function (Q(x) = 1 for x > 0) and 

dQ(k x ,S) J 



dS 




Q(k x ,S) 



(12) 



e(k x ) = cos(k x ) - pk 



(14) 



For convenience we also absorb S into J: J = JS, and drop it 
from the following discussion. We then compute the one loop 
energy as a function of J, fi, and a, and determine the propen- 
sity for magnetic ordering in various directions. From (8) we 
see that the x and z directions become equivalent at one loop, 
so it suffices to set J x = and work with nonzero J y and J z . 
Symmetry considerations show that J = is an extremum of 
the one loop energy. In fact, we empirically see that it is a 
global maximum, and that the energy is unbounded from be- 
low, becoming more negative with increasing \J\. This makes 
sense since the saddle-point Hamiltonian treats S as a clas- 
sical field with no dynamics, and larger 15*1 leads to lower 
electronic energy. Physically, we expect that \S\ ultimately 
saturates. For the purposes of the saddle-point approximation 
we pick an appropriate value for \J\ and evaluate the energy 
difference between ordering in the y and z directions: 

SE(a, n) = E(J, 0, a, fj.) - E(0, J, a, fx) 

~ J 2 (cz(a, /i) — c x (a, /i)) (15) 

In figure 2 we plot 5E/\E\ at a specific value of J, and note 
that it is everywhere negative. Although the percentage dif- 
ference in energies is small, it is robustly negative over a very 
large range of physical values of n, a, and J. We conclude 
that in saddle-point approximation the spins prefer to develop 
magnetic order in the (in-plane) y direction. 

We can perform a similar analysis for a two dimensional 
version of (2), obtaining a similar result: the spins prefer to 
order in-plane in the large-TV limit. 

One very interesting part of our analysis is that we find fer- 
romagnetic order developing at weak coupling. The underly- 
ing reason for this is that the band structure in the presence of 
Rashba spin-orbit coupling has a crossing at k = 0; a small 
magnetic moment opens a gap there. This occurs even at arbi- 
trarily weak coupling, if the chemical potential passes through 
this crossing. This is very similar to the case of other Fermi 
surface instabilities, such as the BCS instability or density- 
wave ordering for nested Fermi surfaces. If the chemical po- 
tential does not pass through the k = crossing, then a small 
minimum coupling must be exceeded, as in the case of small 
detuning away from a nested Fermi surface. This scenario 
stands in stark contrast to the usual case of the Stoner insta- 
bility: ordinarily, ferromagnetism does not open a gap at the 
Fermi surface, and only occurs when the coupling exceeds the 
inverse of the density of states. 
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FIG. 2. -(SE)/\E\, where SE is defined in (15), plotted at J = 0.4 
as a function of chemical potential p and spin-orbit coupling a for a 
one dimensional conducting wire aligned in the y direction sitting on 
the 2 dimensional xy plane. We note that this quantity is everywhere 
negative, implying a propesity towards magnetic ordering in plane, 
perpendicular to the wire i.e. along the x direction. Other values of 
J give similar results. 
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FIG. 3. (color online) Ground state phase diagram of the model 
Hamiltonian in Eq.(16) at filling p — 1/6, determined by accurate 
DMRG simulations with system size up to N = 192 sites. Changing 
coupling parameters J and aa -1 , three different phases are found, 
including the fully polarized phase, the partially polarized phase, as 
well as the easy-plane partially polarized phase. Here J c hain = O.lt 
and U = 40t. 



IV. DMRG SOLUTION OF THE ID LIMIT 

We now consider a ID lattice Hamiltonian describing itin- 
erant electrons coupled to localized impurity spins: 

H = -t^2(cf a c i+ ia + /i.e.) + J^Si- cf a a a pc il3 

— Jchnin^^Si ■ Sj+i (16) 



i 



Here cf (Cj Q ) is the electron creation (annihilation) opera- 
tor with spin index a = (t, ■!•) at s i* e i s tne & = h 
spin operator, representing the localized magnetic moment; t 
denotes the nearest neighbor (NN) tunneling matrix element 
(henceforth we set t — 1) and a is the lattice constant. J 
is the Kondo coupling between localized magnetic moments 
and itinerant electrons, J c hain is the ferromagnetic exchange 
coupling between NN localized magnetic moments, and a is 
the spin-orbit coupling. U is the on-site Hubbard repulsion 
for the itinerant electrons, and we have also included Zeeman 
terms for both the localized spin and itinerant electron spin 
along the z and y directions. In the Hamiltonian (16), impu- 
rity spins have their own dynamics and the ground-state of the 
is determined by taking into account both electron and spin 
degrees of freedom on equal footing. If we use mean-field 
approximation for impurity spins, i.e. (Si) = S(xi), and ne- 
glect electron-electron interaction, we obtain the Hamiltonian 
(2) considered in the previous section. As we show below, our 
conclusions regarding preferred magnetization persist in the 
strongly interacting limit U 3> t. 

In our picture, most of the electrons required by the po- 
larization catastrophe argument become localized d xy spins, 
with a much lower density of itinerant electrons. Therefore, 
we focus on the low density case, taking p — 1/6 for the sake 
of concreteness in most of our calculations. The coupling pa- 
rameters are taken to have values U = 40t and J c ham = t/10. 
(However, our results are not very sensitive to the value of 
•/chain-) We then map out the phase diagram as a function of 
ao _1 and J. By the arguments given in Sec. II, we expect 
that era -1 ~ t/40 and J m — 0.3t. However, given that there 
is some uncertainty in these parameters, it behooves us at this 
stage to see how the physics of our Hamiltonian changes as 
we vary them. We have included small magnetic fields in the 
y- and z— directions. These fields are necessary to break time- 
reversal symmetry and rotational symmetry about the S^-axis; 
otherwise, we would necessarily find (S y ) = (S z ) = 0. In an 
infinite system, but in the absence of these symmetry-breaking 
fields, the system can spontaneously choose to order along 
S y or —S v or it could spontaneously pick a direction in the 
S z — S x plane (if it orders at all). Employing the unbiased 
density matrix renormalization group method, we determine 
the ground-state phase diagram of the system (16). 

We find the following phase diagram, depicted in Fig. 3. For 
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FIG. 4. (color online) Relative polarization M/M s as functions of 
aa' 1 and J, along y direction (M = S v ) and z direction (M = 
S z ), for the system in Eq.(16) at filling p — 1/6 and system size 
N = 192 sites. Relative polarization M/M a for J = 0.3 in (a), 
aa" 1 = 0.0 in (b), aa' 1 = 0.1 in (c), and aa" 1 = 0.2 in (d). Here 
Jchain = O.lt, U = 40t, and M s is the saturated magnetization. 
Note that the polarization is induced by applying a small magnetic 
field h y = 0.005t along y direction or h z — 0.005i along z direction 
separately. 



J large and negative (i.e. a ferromagnetic coupling between 
localized and itinerant spins) and aa' 1 large, the system is 
partially-polarized and the moment points in the y-direction 
(i.e. in plane, but perpendicular to the nanowire, and corre- 
sponding to the light blue region in the upper left corner of 
Fig. 3). For J < but lying below the blue phase boundary 
in Fig. 3, we have a fully-polarized phase in which the spins 
can point equally-well in any direction (the grey region near 
the middle of Fig. 3). For J > 0, there is a partially-polarized 
phase at small J or small aa' 1 (narrow purple region). In this 
region, the spins can point can point equally-well in any di- 
rection. Finally, if J is positive and either J or aa' 1 is large, 
then the system is partially-polarized and the spins point in the 
^-direction (light blue region in upper right of Fig. 3). 

We now examine the magnetically-ordered state in more 
detail. The occupation number n(k) = (clc^ + c|c^) clearly 
shows a filled Fermi sea with approximately one electron in 
each occupied state, with \k\ < kp, as illustrated in Fig. 5a, at 
the parameter values listed in Fig. 5b. The Fermi wavevector 
is consistent with kp = tt/3, in agreement with Luttinger's 
theorem. The structure factor p(q)p(—q) shown in in Fig. 5b 
has cusps at ±2kp, from which we can more precisely extract 
the Fermi wavevector kp — tt/3. The occupation numbers 
and Luttinger volume are consistent with the chemical poten- 
tial depicted in Figure 1 . There are two Fermi points and there 
is a single state at each Fermi point because the spin is locked 
to the momentum. Such a ID electron gas is often called a 
"helical wire". (Note that our system is not simply fully spin- 
polarized; we have checked that, at these parameter values, 
is ever equal to one, so the spins are 
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FIG. 5. (color online) The Fermi surface of our ID model at three 
points in the phase diagram, (a) The occupation number n{k) = 
(clef + cjcj_) , which shows the region of momentum space occupied 
by the filled Fermi sea. (b) The equal-time density-density correla- 
tion function, which has singularities at 2kF- 



not polarized in the z-direction; they are also not equal to each 
other, so the spins are not polarized in the x- or y- directions.) 
Therefore, an odd number of bands (in fact, just one) crosses 
the Fermi surface; as we discuss in Section V, this means that 
the system is primed for the development of topological su- 
perconductivity. 



V. PROXIMITY-INDUCED SUPERCONDUCTIVITY 

We propose the following picture for superconductivity at 
the LAO/STO interface. Our starting point is superconductiv- 
ity in STO, which occurs when insulating STO is doped. We 
assume that insulating STO has small islands or droplets of lo- 
cal superconductivity, which are too far apart and too weakly- 
coupled to develop long-ranged superconducting order. We 
suppose that these droplets are caused by unintentional local 
defects in STO. When STO is doped, the islands of local su- 
perconductivity grow in size and become more strongly cou- 
pled, until long-ranged superconducting order sets in. How- 
ever, the presence of an interface with LAO changes matters. 
Itinerant electrons at the interface can mediate a coupling be- 
tween superconducting droplets in STO that are close to the 
interface. As we will show, this can enable superconductivity 
to develop even when the droplets are too weakly-coupled to 
percolate across STO on their own. 

We show how this can occur with a calculation in a simpli- 
fied model. We suppose that there are some superconducting 
droplets in STO that are near the LAO/STO interface. In each 
droplet, a single-particle gap is assumed to be well-developed, 
but the interactions between the droplets are assumed to be too 
weak for superconducting order to set in. AID channel at the 
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interface couples to a subset of these droplets, which form a 
linear array. The ID channel induces interactions between the 
droplets, so that the linear array can be modeled as a ID spin- 
gapped electron system, which we assume to be just slightly 
on the disordered side of the Kosterlitz-Thouless transition. 
(Since doped STO superconducts, this is a reasonable approx- 
imation.) The coupling between this system and a ID channel 
at the LAO/STO interface can nudge the system into the basin 
of attraction of the quasi-long-range ordered superconducting 
phase on the other side of the Kosterlitz-Thouless transition. 
We thereby see how proximity to a metallic interface can sta- 
bilize long-ranged or quasi-long-ranged superconducting or- 
der. 

Once superconductivity is established in STO, it is induced 
at the interface by the proximity effect. This mechanism of es- 
tablishing superconducting order works for generic conduct- 
ing ID channels. However, the interesting scenario for us 
occurs when, due to Rashba spin-orbit coupling and ferro- 
magnetism, the ID channel realizes the helical band struc- 
ture in Fig 1, and the induced superconductivity is topolog- 
ical. Hence, we begin with a bosonized helical ID channel at 
the interface, coupled to an array of isolated, i.e. effectively 
zero dimensional, droplets: 

S=^K W J dxdr [{d x 9) 2 +v F - 1 (d T 9) 2 ] 

+ ^E(^) 2 + A ^E {^ 9{X]) e ~^ 3 + c - c -) 

3 j 3 

(17) 

Here, vp is the Fermi velocity in the ID channel and K w is 
its Luttinger parameter. The ID channel is assumed to have 
repulsive interactions, so K w < 1. The factor of V2 in the ex- 
ponent in the third line has been inserted for later convenience 
so that it agrees with the convention for the charge boson of a 
spin-gapped electron system (see, for instance, Ref. 30). The 
droplets are assumed to have an average spacing a, an average 
Josephson coupling Ap to the ID channel, and an average 
charging energy U. We will neglect random variations in the 
spacing between droplets, in the Josephson couplings, and in 
the charging energies, and simply set Xj = ja, Apj = Ap, 
and Uj = U . Random variations in these parameters are cer- 
tainly important in the physical system but are an unnecessary 
complication for this calculation. Note that if Ap = 0, then 
the droplets are completely decoupled from each other and 
there is no quasi-long-ranged superconducting order. 

We now integrate out fluctuations of 6 at length scales 
shorter than I. This generates a coupling between droplets. At 
length scales much longer than I, we can take the continuum 
limit for the array of droplets, thereby leading to the following 
effective action: 

S = 2tt j dxdT { VpKw ^ dx0)2 + VF ~*(W\ 

+ vk p [(d x e p ) 2 + v- 2 (d T e p ) 2 ] 

- — cos (26- V29 p ) ) (18) 



The droplets are now effectively described by a ID wire with 
a spin gap; 6 P is the charge boson for such a wire. Since we 
assume that there is, initially, a weak interaction between the 
droplets and the ID channel, we assume that K p < 1. In other 
words, although the wire has a spin gap, we do not assume 
that it can superconduct without further mediation on the part 
of the ID channel at the interface. Indeed, the parameters K p 
and v can be related to the effective superfluid stiffness p s 
and compressibility k of the array: K p = 2-k^J 'A w p s K and 
v = \J ' A w p s / ' k with being A w the cross-sectional area. We 
assume here that the superfluid stiffness is such that K p < 1. 

If the two velocities were equal, v — vp, one could analyze 
the model by forming the combinations 0± = {y2Q p ± 20) /2. 
In terms of new variables 6± the action reads 

S=^Jdxdr ((^ + ^) [(d p 6 + ) 2 + (d p 6_) 2 ] 

- 2 («f - ^) duB+dJ- - 2y cos (20_)J . (19) 

Here, we have rescaled the time coordinate by v and have in- 
troduced the dimensionless parameter y = Apa/2v. We can 
now see, at a heuristic level, how the coupling between the 
droplets and the ID channel can stabilize quasi-long-ranged 
superconducting order. Let us suppose, for a moment, that 
the coupling y is relevant. Then is pinned, and we can ig- 
nore its fluctuations. Then we are left with 8 + , which exhibits 
algebraically-decaying superconducting order. This order is 
stable if weak impurity-backscattering or, equivalently, vor- 
tex tunneling is irrelevant. Because the wire is helical, we can 
only tunnel — vortices 30 , i.e. 9 —> 9 + 2n. Since 6- is pinned, 
this means that 9 p must also wind 9 p — >• 9 p + 2tt\J2. Conse- 
quently, in such a process, 9 + — s- 9 + + 4tt. The operator that 
accomplishes this is sin(4<^+), where [<j>+(x),d x 6+(y)] = 
i / w8{x — y). This operator is irrelevant if 2K p + K w > 1, 
and when this inequality is satisfied, the system exhibits quasi- 
long-ranged superconducting order. Note that this can be sat- 
isfied even if K p < 1 and K w < 1. So two systems, neither of 
which could sustain superconductivity on their own, can de- 
velop superconductivity when in proximity to each other. The 
key to this is the topological nature of the superconductivity: 
since only vortices can tunnel through a helical wire, the 
stability condition is less strict than for an ordinary supercon- 
ductor 41 . 

To support the aforementioned scenario, we need to show 
that that Cooper-pair tunneling term is relevant. The corre- 
sponding RG equation for y is: 

^ = [2-( K >)- 1 ]y, (20) 

where (K 1 )- 1 = ^ + If (K 1 )- 1 < 2, y will grow 
from the initial small value y(0) <C 1 to y(l) ~ 1 at which 
point 6>_ gets pinned. Given that K' also flows under RG, we 
need to compute its flow equation and complete the system of 
RG equations for this model. To do that we rewrite Eq. (19) 
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in the following form: 



S = 



(W 5 



2K + -{d fi + d IM 6-) - 2ycos(26_) (21) 



where, initially, K+ = K- = 

K, 



K, = 



Kj, 

4 1 2 ' "H — 4 

The reason that we have introduced three couplings 
K + , K- , when there are, seemingly, only two cou- 
plings K w and K p is that the RG flow for this theory will 
carry the system away from the initial point K + = A'_. A 
real-space RG calculation yields the following equations for 



dK + dK+. 



= 0. 



(22) 



The coupling K' can be expressed in terms of K_, K + and 
A'_| (or equivalently in terms of K- and the initial parame- 
ters K p and K w since and A + do not flow): 



K' = K_ — K%_/Kj i 



A_ 



{2K p - K w ) 2 
A(2K p + K w ) • 



(23) 



Since A' monotonically depends A_, the growth of A_ un- 
der the RG flow results in 



dK' 



V 



(24) 



Thus, the j/-coupling becomes more and more relevant, and 
eventually pins 6>_ as assumed above. 

We now consider more general case of unequal velocities 
v 7^ v F . Proceeding as before, we find that Eq. (19) can be 
written as 



We now sketch the real-space RG procedure. We start out 
by integrating out short-distance modes, but allowing arbitrar- 
ily short times. Thus, we have an effective action in which 
there are modes 6(k,ui) with \k\ < a -1 and < oo, where 
a -1 is the momentum cutoff. Following a real-space RG ap- 
proach , we integrate shells a < r < as while keeping 
time integrals unconstrained and eventually rescale r — > sr, 
t st. Here s = e dl . Our RG procedure involves calculat- 
ing the correlation function ( e 2i M Xl ' T1 ) . e -^0-{x 2 ,r 2 )y To 
do it safely one has to normal order the exponent: 

e 2»Ml) . e -2«M2) _. e 2i[0-(l)-0_(2)] . e -2([0_(l)-e_(2)] 2 ) 

(30) 

where 6(1) = 6(x\,tx), and the average (...) is computed 
with respect to the bare action (y = 0) defined in Eq. (25). 
Initially, when A = K + and v — v + , the correlation function 
in Eq. 30 can be easily calculated 



-2([Mx,t)-M0,0)] 2 ) 



-J— + -L- 

a 2K p ' K w 



(x 2 + V 2 T 2 ) 2K e (x 2 + V 2 F T 2 ) 



(31) 



However, in general it is a complicated function of the cou- 
pling constants (26) as well as A and v . One can show that 
at the tree-level the RG equation for y becomes 



dt 



= [2-(K'y l ]y 



(32) 



where, again, (A') 1 

We now compute the RG equations for this model at one- 
loop level. Using Eqs.(30) and (31), we find 



5 = — / tl.r (It 
2tt 



(R ^(d^^+V- 1 ^) 2 ) 

+ K+ (i+(d x 6+) 2 + v+ 1 (d T 6 + ) 2 ) 
2K+- (v+-(d x 6+)(d x 6-) +v^_(d T 6 + )(d T 6-)) 



2j/cos(20_) (25) 



where the coupling constants are defined as 
1 



A, = 



(K w v F + 2K p v) (K w v + 2K p v F ) (26) 



v+ 



2K p v + K w v F 



' VVp 



2K p v F + K w v 



(27) 



v+- 



1 



(K w v F - 2K p v) (K w v - 2K p v F ) (28) 



' VVp 



K w v F — 2K p v 
K w v — 2K p v F 



(29) 



The parameters A and v are initially equal to K + and v + , 
respectively, but they flow under RG as explained above. 



'VVp 



(33) 
(34) 



where y = Apa/y/vvp and the dimensionless function 



U(k) 



1 

8^ 



dz- 



(k- 1 + z 2 ) 2K p (k + z 2 ) 



(35) 



As follows from Eqs.(33), the RG equations for A and v are 
given by 



dK 


_v 2 


H 


2 


dv 


1 


d£ = 


2K 



L/„p) + ^/ 2 p) 

vvp \ v J v \ v I 



\JVVp 
IVVpfa 



/VVp 



(36) 
(37) 



Thus, according to these Kosterlitz-Thouless-type RG equa- 
tions, Eqs. (32),(36) and 37, we see that K grows. The pa- 
rameter K' has a complicated dependence on K and v which 
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VI. DISCUSSION 




FIG. 6. Dependence of K' on the flow parameters K and v. Here 
we used K m = 0.8, K p — 0.2 and v — 0.5vf- The function K' is a 
monotonically increasing function of K. 



follows from Eq. (31). However, at small initial velocity mis- 
match \8v\ =\v — vp\ <C vf, one finds 



(K>) 



K 



(2K P - K w f 



A(2K p + K w )_ 
32K(K 2 - 4K 2 )K p (3K n 



(38) 



2K P ) Sv 2 



A(2K p + K w )K-(2K p -K w ) 2 



One can see that the 5v 2 -correction is quickly decaying with 
K and thus do not change qualitatively our results obtained for 
the v = vp case, cf. with Eq. (23). In general, we find that K' 
is a monotonically increasing function of K, see Fig. 6. Thus, 
the growth of K implies the growth of K' . Therefore, once 
relevant, y will grow to strong coupling y(l*) ~ 1 and pin 
0_. In this case, 0_ drops out from Eq. (25) and the effective 
action now reads: 

s= *h j " dxdT ^(n(s(O(dx0+) 2 +6(n~ 1 (dT0+) : 

(39) 

Following the argument that we used for v = Vp, we observe 
that quasi-long-ranged order is stable so long as flux hc/e vor- 
tex tunneling is irrelevant, i.e. when sin(4</> + ) is irrelevant. 
This occurs when AK > 1. Note that this can be satisfied 
even if K p < 1 and K w < 1. As in the equal velocity case, 
two systems, neither of which could sustain superconductivity 
on their own, can develop superconduct when in proximity as 
a result of the helical nature of one of the systems. 

We thereby arrive at the model of Ref. 30: a ID chan- 
nel that is proximity-coupled to a quasi-long-range-order su- 
perconducting wire. As shown there, such a wire supports 
Majorana zero modes. We also expect our results discussed 
in this section to apply to multichannel nanowires with an 
odd number of occupied subbands coupled to superconduct- 
ing droplets, see, e.g., Refs. A2-A4. 



In this paper, we have adopted the point of view that 
SrTiC>3 has the seeds of both magnetism and superconduc- 
tivity. However, these local tendencies only come to fruition 
when brought into contact with a metallic layer or ID channel. 
We have focussed on the latter case, for reasons of tractabil- 
ity as well as potential relevance to the experiments of Refs. 

and , but we believe that our general mechanism works 
in 2D as well. We have shown that local moments in SrTiO,3 
that are near the LAO/STO interface can order ferromagneti- 
cally, as a result of their interaction with mobile electrons at 
the interface. We have also shown that droplets of local su- 
perconductivity in STO - which would interact too weakly to 
develop superconducting order if left to their own devices - 
can develop superconducting order as a result of their inter- 
action with mobile electrons at the interface. Finally, we have 
noted that the interface electrons can form a topological super- 
conducting state as a result of their proximity to ferromagnetic 
and superconducting order. 

We have shown that our model leads to ferromagnetism by 
two different calculations: a large- A calculation and a DMRG 
calculation. Both calculations find a ferromagnetic state with 
spins pointing in the plane, i.e. along y-axis with x being the 
direction along the wire. (In the DMRG calculation, the polar- 
ization may be either partial or full, depending on the strength 
of the spin-orbit coupling.) Interestingly, our large- N calcula- 
tion finds a ferromagnetic state even at weak coupling, which 
is a feature of the band structure in the presence of Rashba 
spin-orbit coupling. AID wire with sufficiently strong spin- 
orbit coupling and Zeeman field will form a helical wire. Our 
calculations - both large- A*" and DMRG - show that our model 
gives rise to a helical wire. Remarkably, recent transport mea- 
surements can be interpreted as evidence that ID channels at 
the LAO/STO interface are helical wires . 

In a helical wire, it is possible for s-wave superconductiv- 
ity to coexist with a magnetic moment. As a result our model 
allows for a proximity coupling between s-wave supercon- 
ducting droplets in STO and ferromagnetic electrons at the 
LAO/STO interface. We have analyzed our model for super- 
conductivity by mapping it to the theory of a single ID bo- 
son in the vicinity of the Kosterlitz-Thouless transition. We 
find that such a model could be on the disordered side of 
the Kosterlitz-Thouless transition for very weak coupling be- 
tween the droplets and a ID channel at the interface but it 
could be on the ordered side of the transition if the coupling is 
sufficiently strong. 

Our results on magnetism and superconductivity imply that 
the superconducting state of a ID channel at the LAO/STO in- 
terface is in a topological superconducting state. This actually 
stabilizes the system against quantum phase slips: 2tt phase 
slips are forbidden, and only Air phase slips, which are less 
relevant in the RG sense, could disrupt the superconductivity. 
Furthermore, this topological superconducting state supports 
Majorana fermion zero modes, whose presence would lead to 
a 47r-periodic ac Josephson effect 24,29 . 

The main goal of our work has been to show why mag- 
netism occurs, why superconductivity occurs, and why they 
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FIG. 7. A schematic phase diagram for the LAO/STO interface as a 
function of carrier concentration. 



portant target for further investigation. 

We note that our calculations rely heavily on simplifying 
features of one dimensional systems - the applicability of 
DMRG calculations to the magnetic ordering of the system 
and the applicability of bosonization to the superconducting 
ordering of the system. It would be interesting to give a fully 
two-dimensional analysis of a model similar to ours. Further- 
more, it would be useful to investigate other related materials 
exhibiting the same phenomena. In particular, recent work 
on epitaxially grown GdTi03-SrTi03 interfaces indicates 
that ferromagnetism and superconductivity can also coexist 
in such systems. Being much cleaner than LAO/STO from 
a materials point of view, these interfaces might provide an 
attractive environment for the investigation of the ideas pro- 
posed in this paper. 



can coexist. However, magnetism is probably not found at all 
carrier concentrations above the metal-insulator transition . 
Neither is superconductivity. A plausible schematic phase di- 
agram is given in Fig. 7, based on the phase diagram in Refs. 

and . Therefore, it is also important to understand when 
and why they do not occur. If the magnetic order is too strong, 
so that the spins are fully polarized, then our mechanism does 
not work. This may explain why superconductivity is sup- 
pressed at large carrier concentration. However, a more de- 
tailed understanding of the phase diagram is definitely an im- 
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